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Measurements have been carried out over the temperature range 380°-525°K using a four-cell hot-wire 
I thermal-conductivity apparatus. At the lower temperatures the thermal conductivities of H 2 O and DjO are / / j 

I identical, but at the highest temperature the conductivity of D 2 O is larger by 8-9 parts per thousand. Con- \~sCo( 

/ ductivities of the equimolar mixtures lie midway between the conductivities of the pure components. Results 
\ are analyzed to test a recent hypothesis that a resonant exchange of rotational quanta causes the thermal 
conductivity of a polar gas to seem anomalously low in relation to its viscosity. It is concluded that the inter- 
change of rotational energy between neighboring dipoles may be important, but that exact resonant is 
‘ ~ 7 not required, at least in the case of water vapor. 



INTRODUCTION 


O UR understanding of heat conduction in poly- 
atomic gases has been considerably enhanced by 
the recent theoretical studies of Mason and Monchick. 1 
Starting with the formal kinetic theory of Wang Chang 
and Uhlenbeck 2 and Taxman, 3 they have derived ex- 
plicit expressions for the thermal conductivity of poly- 
atomic gases. By systematically including terms in- 
volving inelastic collisions, they derived the modified 
Eucken expression*"" 6 as a first approximation and, as 
a second approximation, an expression dependent on 
the relaxation times for the various internal degrees 
of freedom. Their result may be written 

<« 

Here A is the thermal conductivity, 17 is the viscosity, 
M is the molecular weight, p is the density, and D m t is 
the average coefficient for the diffusion of internal 
energy, while C* trang and C- mt are the translational and 
the internal contributions to the molar heat capacity. 
Z k is the number of collisions for relaxation of the £th 
internal mode, while C* is the heat capacity associated 
with that mode. 

The first two terms in Eq. (1) are simply the modi- 
fied Eucken approximation, while the third term is 
important only for small collision numbers. In small 
rigid polyatomic molecules these are associated with 
rotational relaxation. (In flexible molecules vibrational 
relaxation times may also be small.) 

For nonpolar gases Dmt is commonly assumed to be 
equal to the self-diffusion coefficient. With this assump- 




tion O’Neal and Brokaw 7-9 have analyzed experimental 
recovery-factor measurements (in essence a direct de- 
termination of the quantity \M/i)C v ) so as to obtain 
collision numbers for rotational relaxation. Since colli- 
sion numbers for nitrogen, oxygen, carbon dioxide, 
methane, and hydrogen were found generally in accord 
with values obtained by other experimental methods, 
it was concluded that the approximate theory embodied 
in Eq. (1) is substantially correct. Furthermore, it 
seems that these collision numbers are determined by 
the following factors: (1) the molecular mass distribu- 
tion, (2) the strength of the intermolecular attractive 
forces, and (3) the molecular asymmetry. 

The thermal conductivities of highly polar gases 
such as hydrogen fluoride, water, and ammonia ap- 
pear to be anomolously low in relation to their viscosi- 
ties. Mason and Monchick 1 suggest that this effect is 
largely due to a resonant exchange of rotational energy, 
presumed probable on grazing self-collisions of polar 
molecules. Hence, a grazing collision with exchange is 
equivalent to a head-on collision without exchange in- 
sofar as the transport of the rotational quantum is 
concerned. Thus, the diffusion coefficient for internal 
energy, D mt in Eq. (1), is smaller than the self -diffusion 
coefficient, Du, and is given by the expression 


Dint= = Dn/(l+3), 


( 2 ) 


where 8 is a correction term calculated from the theory 
of resonant collisions. 

Mason and Monchick 1 give expressions for 8 for 
linear dipoles and also for several types of symmetric 
tops. The salient factors determining 8 are apparent 
in their result for spherical tops: 


1 E. A. Mason and L. Monchick, J. Chem. Phys. 36, 1622 
(1962). 

* C. S. Wang-Chang and G. E. Uhlenbeck, University of 
Michigan Engineering Research Rept. No. CM-681 (July 1951) . 

* N. Taxman, Phys. Rev. 1 10, 1235 ( 1958) . 

4 S. Chapman and T. G. Cowling, The Mathematical Theory of 
Non-Uniform Gases (Cambridge University Press, New York, 
1952), p.238. 

*K. Schafer, Z. Physik. Chem. B53, 149 (1943). 

6 J. O. Hirschfelder, J. Chem. Phys. 26, 282 (1957). 


a _5r (0l) rAjt? 

16 k 'Y\(><T D kT\ 




(3) 


L(Cint/R)/ } J' 

Here (a 4 ) is the mean value of a dimensionless quantity 


7 C. O’Neal, Jr., and R. S. Brokaw, Phys. Fluids 5, 567 (1962). 

8 R. S. Brokaw and C. O’Neal, Jr., Symp. Combust. 9th, Cornell 
Oniv., 725 (1963). 

» C. O’Neal, Jr., and R. S. Brokaw, Phys. Fluids 6, 1675 (1963). 
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involving the rotational quantum numbers, taken to 
be 0.44. Furthermore, m is the dipole moment, h is 
Planck’s constant, an is the kinetic-theory diameter 
characteristic of diffusion in the notation 

of Ref. 10), k is the Boltzmann constant, m is the 
molecular mass, R is the universal gas constant, and 
I is the molecular moment of inertia. [It might be 
noted that <5 for linear dipoles is given by an expression 
identical to Eq. (3) but for the factor 57r/16.] 

As already noted, the collision number is dependent 
on the molecular mass distribution. This dependence 
may be characterized as 

Z t «r'*4I/m<r*, (4) 

where <r v is the viscosity collision diameter ( = o-[S2 (2>2) *J 
in the notation of Ref. 10) . 

It occurred to us that Mason and Monchick’s postu- 
late regarding the importance of the resonant exchange 
of rotational quanta might be tested by measuring the 
thermal conductivities of deuterated and ordinary polar 
gases; H 2 0 and D 2 0 were chosen because they were 
readily available. 

From Eqs. (l)-(4) it is seen that isotopic substitu- 
tion of deuterium for ordinary hydrogen will affect 
the thermal conductivity through increases in molecu- 
lar weight (10%), internal heat capacity (10%-15%) , 
and moments of inertia (80%-100%). Thus the mo- 
ments of inertia are profoundly altered. In the absence 
of the resonant phenomenon we would expect the ther- 
mal conductivity of H 2 0 to be somewhat larger than 
that of D 2 0. This is a consequence of the smaller mass of 
ordinary water [see Eq. (1) — viscosities are only pro- 
portional to M*] and also of the smaller moments of 
inertia which lead to a larger collision number for 
rotational relaxation [see Eq. (4)]. These factors 
should outweigh the effect of the somewhat greater 
heat capacity of D 2 0. 

On the other hand, the resonant correction depends 
inversely on moment of inertia [Eq. (3)] and, hence, 
is larger for H 2 0. As a consequence, if the resonant 
exchange of rotational quanta does indeed occur, the 
conductivities of H 2 0 and D 2 0 should be very nearly 
equal; the conductivity of D 2 0 may actually be slightly 
larger. 

Mason and Monchick’s postulate regarding resonant 
exchange of rotational energy can be subjected to a 
further qualitative test. The effect is based on exact 
resonance; that is, an exact matching of rotational 
energy levels. Consider now an equimolar mixture of 
H 2 0 and D 2 0. In such a gas one-half of the collisions 
are self-collisions of H 2 0 or D 2 0 in which resonant 
exchange of rotational quanta may occur. The remain- 
ing collisions are between H 2 0 and D 2 0. In these col- 
lisions there is no matching of rotational levels so there 
should be no resonant exchange. Consequently, we 

10 J. 0. Hirschfelder, C. F. Curtiss, and R. B. Bird, Molecular 
Theory of Gases and Liquids (John Wiley & Sons, Inc., New York, 
1954). 
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might expect such a mixture to have a thermal con- 
ductivity larger than either of the pure constituents. 
This argument must remain qualitative, however, be- 
cause we have as yet no analogous theory for mixtures. 

The present experiments were initiated to test these 
theoretical notions regarding the thermal conductivity 
of polar gases. (A preliminary result using a rather 
crude apparatus has already been reported. 11 ) Con- 
ductivities of H 2 0 and D 2 0 were measured at 381.2°, 
426.1°, 478.0°, and 525.6°K, while the thermal con- 
ductivity of the equimolar mixture was measured at 
the two highest temperatures only. After completion 
of the present experimentation and data analysis we 
learned of the recent work of Vargaftik and Zaitseva 12 
who have measured the thermal conductivities of H 2 0 
and D 2 0 from 100°-500°C. The present measurements 
are in qualitative, but not quantitative, agreement 
with the Russian work. 

EXPERIMENTAL 

Apparatus and Procedure 

The hot-wire thermal-conductivity apparatus has 
been described in some detail by Coffin and O’Neal. 13 
Briefly, it consists of two pairs of hot-wire cells con- 
nected as elements of a constant-current Wheatstone 
bridge. The voltage unbalance produced when gases 
of unequal thermal conductivity are introduced into 
the pairs of cells is measured. The instrument responds 
to the reciprocal of the thermal conductivity, 

(5) 

Here E is the voltage unbalance with a gas of unknown 
conductivity in one pair of cells and a reference gas in 
the other pair; E Te{ is the voltage with the reference 
gas in both cells. X and X re f are the thermal conductivi- 
ties of the unknown and reference gases, respectively, 
while b is a constant characteristic of the apparatus 
(it may be slightly temperature dependent). For this 
work helium and argon were used as calibrating gases. 
The conductivities assumed for helium and argon are 
discussed in the next section. 

H 2 0 and D 2 0 vapors were supplied to the cells at a 
pressure slightly above atmospheric from ampoules ther- 
mostated in a bath at 104°C (a boiling salt solution) . 
Connecting lines were heated with electrical heating 
tape to prevent condensation. The gas handling system 
was in other respects conventional. 

In addition to the measurements on ordinary and 
heavy water vapors, the thermal conductivity of nitro- 
gen was determined at all four temperatures. Also, air, 
oxygen, and carbon dioxide were measured at 478.0°K. 

11 Reference 8, p. 721. 

12 N. B. Vargaftik and L. S. Zaitseva, Inzh.-Fiz. Zh., Akad. 
Nauk Belorussk. SSR 6, (5), 3 (1963). 

13 K. P. Coffin and C. O’Neal, Natl. Advisory Comm. Aeronaut. 
Tech. Note 4209, Washington, D.C., 1958. 
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Table I. Experimental emf measurements 
relative to argon (mV). 


Temperature, °K 


Gas 

381.2 

426.1 

478.0 

525.6 

He 

5.030 

5.167 

5.186 

5.208 

n 2 

1.730 

1.800 

1.791 

1.802 

d 2 o 

0.806o 

1.126 6 

1.413s 

1.640s 

h 2 o 

0.808* 

1 . 128s 

1.386o 

1.604* 

0.5H 2 0-0.5D 2 0 



1.400o 

1.626* 

O* 



2.091 


co* 



0.980 


Air 



1.861 
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In general, nitrogen served as the reference gas except 
that D 2 0 and the equimolar H 2 0-D 2 0 mixture were 
measured against ordinary water. Thus these small 
emf’s were measured with enhanced precision. 

Experimental voltages relative to argon are recorded 
in Table I. An extra significant figure is given for the 
isotopic waters to reflect the measured voltage differ- 
ences among them. The voltages relative to nitrogen 
were reproducible to about 0.001 mV. 


Fig. 1. Departure of experimental thermal conductivities from 
tabulated values for helium. (Solid symbols are conductivity 
measurements; open symbols calculated from viscosity data.) 

est this quotient is greater than unity by 1-3 parts 
per thousand for argon and 3-4 parts per thousand 
for helium. 

In essence, large scale plots of experimental thermal 
conductivity were constructed and smooth curves 
faired through the data, 16 ^ 19 giving the greatest weight 


Thermal Conductivities of the Calibrating Gases 
Helium and Argon 


Thermal conductivities for the calibrating gases were 
assigned after considering experimental viscosity and 
thermal conductivity data on helium and argon. Be- 
cause helium and argon are monatomic gases, it is 
possible to compute conductivity from viscosity through 
the rigorous theoretical relation 14 


15 R 

4 M 71 /„ (3) ' 


( 6 ) 


The functions /x (3) and /, (3) are very slowly varying 
functions of temperature that differ only slightly from 
unity. Equation (6) is rigorous through the third 
Chapman-Enskog approximation and has been experi- 
mentally verified for helium and argon within a few 
parts per thousand. 7 The factor /x (3) //i» (3) was calculated 
assuming an exponential-6 force law with appropriate 
force constants 16 ; over the temperature range of inter- 


Table II. Thermal conductivities of helium and argon assumed 
for calibration purposes (peal cm -1 sec -1 °K“ 1 ) . 


Temperature, 

°K 

Helium 

Argon 

300 

370.9 

42.42 

381.2 

437.4 

51.26 

400 

452.3 

53.27 

426.1 

472.8 

55.82 

478.0 

512.8 

60.79 

500 

529.5 

62.85 

525.6 

548.6 

65.12 

600 

603.4 

71.51 


14 Reference 10, Eq. (8.2-32), p. 534. 

15 Reference 10, pp. 181, 1173-4. 
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Fig. 2. Departure of experimental thermal conductivities from 
tabulated values for argon. (Solid symbols are conductivity 
measurements; open symbols calculated from viscosity data.) 


16 Helium viscosity: J. Kestin and W. Leidenfrost, Physica 25, 
537, 1033 (1959) ; R. Wobser and F. Muller, Kolloid-Beih. 52, 165 
(1941) ; M. Trautz and H. E. Binkele, Ann. Physik 5, 561 (1930) ; 
M. Trautz and I. Husseini, Ann. Physik 20, 121 (1934). 

17 Helium thermal conductivity: L. S. Zaitseva, Zh. Tekhn. Fiz., 
29, 497 (1959) [English transl.: Soviet Phys. — Tech. Phys. 4, 444 
(1959)]; W. G. Kannuluik and E. H. Carman, Proc. Phys. Soc. 
(London) B65, 701 (1952); H. L. Johnston and E. R. Grilly, J. 
Chem. Phys. 14, 233 (1946). 

18 Argon viscosity: Kestin and Leidenfrost, Trautz and Binkele, 
Wobser and Miiller, see Ref. 16; A. G. DeRocco and J. O. Hal- 
ford, J. Chem. Phys. 28, 1152 (1958); C. F. Bonilla, S. J. Wang, 
and H. Weiner, Trans. ASME 78, 1285 (1956) ; T. Maluta, Rev. 
Phys. Chem. Japan 27, 16 (1957). 

19 Argon thermal conductivity: Zaitseva, Kannuluik, and Car- 
man, see Ref. 17; K. Schafer, Dechema Monograph 32, 61 
(1959); W. C. Gardiner and K. Schafer, Z. Elektrochem. 60, 588 
(1956); W. F. Shottky, Z. Electrochem. 56, 889 (1952); L. A. 
Bennett and R. G. Vines, J. Chem. Phys. 23, 1587 (1955); R. G. 
Vines, Australian J. Chem. 6, 1 (1953); A. Michels, A. Botzen, 
A. S. Friedman, and J. V. Sengers, Physica 26, 121 (1956); F. G. 
Keyes, Trans. ASME 76, 809 (1954) . 
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Table III. Experimental thermal conductivities 
(/ucal cm' 1 sec -1 °K“ l ) . 


Gas 

381.2 

Temperature, °K 
426.1 478.0 

525.6 

D 2 0 

59. 7 2 

69. U 

80. 0i 

90. 1 6 

h 2 o 

59 . 7 4 

69. 1 4 

79.52 

89. 3* 

0.5H 2 0-0.5D 2 0 



79. 7 7 

89. 8 6 

n 2 

73.6 

80.6 

87.4 

93.7 

0 2 



94.3 


co 2 



72.9 


Air 



88.9 


a D 2 0 

57.6 

66.8 

78.8 

90.6 

a H 2 0 

59.2 

68.2 

79.1 

90.0 


a Reference 12. 


to values computed from recent viscosity measure- 
ments, since these data have the highest precision and 
accuracy. 

The conductivities selected for calibration, as well 
as values at even increments of 100°K are given in 
Table II. Figures 1 and 2 show the departures of 
experimental thermal conductivities from the tabu- 
lated values for helium and argon, respectively. The 
data indicated as solid symbols represent actual ther- 
mal conductivity measurements, while the open sym- 
bols are derived from viscosity measurements and show 
considerably less scatter. 

Experimental Results 

The thermal conductivities derived from Tables I 
and II by means of Eq. (5) are presented in Table III. 
Again an extra significant figure has been retained for 
the isotopic waters so as to indicate relationships among 
them. The values of Table III are generally in accord 
with previous determinations in the literature. This is 
illustrated in Fig. 3, where the present data on water 
vapor and nitrogen are compared with the measure- 
ments of Vargaftik and Oleschuk 20 and Geib and 



Fig. 3. Comparison of nitrogen and water thermal conductivities 
with literature values. 


20 N. B. Vargaftik and 0. N. Oleschuk, Izv. Vses. Teplotekh. 
Inst. 15 (6), 7 (1946). 


Schafer. 21 The agreement is generally satisfactory, al- 
though our datum for water at the highest temperature 
seems to be a bit low. 

Also included in Table III are thermal conductivities 
of H 2 0 and D 2 0 interpolated from the recent data of 
Vargaftik and Zaitseva. 12 The agreement in the case 
of H 2 0 is entirely satisfactory with deviation from 0.5% 
to 1.3%. In the case of D 2 0, however, there seem to 
be more serious discrepancies with deviations of as 
much as 3.5% at the lowest temperature (deviations 
are systematic rather than random). The apparatus 
used by Vargaftik and Zaitseva 12 had a single hot-wire 
thermal-conductivity cell so that H 2 0 and D 2 0 were 
of necessity measured in separate runs. In contrast, 
in our experiments we have compared H 2 0 and D 2 0 
directly against one another. Consequently, we believe 
our results on D 2 0 are to be preferred — especially with 
regard to the values relative to H 2 0. 


Table IV. Molecular properties of H 2 0 and D 2 0. 



H 2 0 

D 2 0 

“Moments of Inertia, 

1.0224 

1.833 

g cm 2 Xl0 40 

b l .9180 

b 3.841 

2.9404 

5.674 

Molecular weight 

18.0160 

20.029 

c Dipole moment, 

1.844 

1.861 

esu-cmXlO 18 




a A. S. Friedman and L. Haar, J. Chem. Phys. 22, 2051 (1954). 
b Figure axis. 

c L. G. Groves and S. Sugden, J. Chera. Soc. 1935, 971. 


ANALYSIS AND DISCUSSION 

The experimental data on the pure gases were ana- 
lyzed as follows: Eq. (1) was assumed correct and 
used to calculate the collision number Z rot and reso- 
nant correction 8 required to reproduce the experimen- 
tal data. Although we have made no assumption as 
to the absolute magnitudes of these parameters, we 
have assumed that the relative magnitudes may be 
calculated from theoretical expressions in terms of mo- 
lecular masses, moments of inertia, dipole moments, 
etc., of the two molecules (values used are shown in 
Table IV). 

Certain temperature-dependent quantities that are 
needed, such as viscosity and heat capacity, are shown 
in Table V. Viscosities of heavy water were computed 
assuming the viscosity collision cross section of D 2 0 
is 3% larger than that of H 2 0. This was deduced from 
the viscosity measurements of Bonilla, Wang, and 
Weiner 18 and is qualitatively in accord with the some- 
what higher dipole moment of D 2 0. 

As was mentioned in the Introduction, Mason and 
Monchick 1 developed expressions for the resonant cor- 
rection for linear dipoles and symmetric tops. They 


21 H. Geier and K. Schafer, Allgem. Warmetech. 10, 70 (1961). 
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Table V. Data analysis for pure HjO and DjO. 



381.2 

426.1 

Temperature, °K 
478.0 

525.6 

Viscosity of H 2 0, /x P 

•127.4 

b 143.2 

M64.0 

b 184.9 

*(CWR) H 2 0 

1.60069 

1.64771 

1.70781 

1.76743 

d 2 o 

1.75277 

1.83355 

1.93036 

2.02233 

VVi?(==6/504*» 

1.316 

1.320 

1.325 

1.327 

Zrot H*0 

2.7 

2.5 

2.3 

2.0 

D*0 

1.6 

1.5 

1.3s 

1.2 

Shso exptl 

0.48 

0.40 

0.38 

0.36 

calc 

0.64 

0.53 

0.44 

0.38 

Sd*) exptl 

0.16 

0.13 

0.12 

0.11 

calc 

0.21 

0.17 

0.14 

0.12 


* J. Hilsenratb, C. W. Beckett, W. S. Benedict, L. Faro, H. J. Hoge, J. F. Masi, R. L. Nuttall, Y. S. Touloukian, and H. W. Woolley, Tables oj Thermal Proper- 
ties of Gases, National Bureau of Standards Circular 564 (1955). 
b J. Kestin and H. E. Wang, Physica 26, 575 (1960). 
c See Ref. a, Table IV. 
d Reference 23. 


suggest that slightly asymmetric tops can be treated 
in the following manner. If I a is the moment of inertia 
about the dipole axis (in this case, the figure axis of 
the molecule) and Ib and Ic are the other two mo- 
ments, one uses the symmetric top formulas replacing 
Ib by ( IbIc )*. From Table IV it is clear that H 2 O 
and D 2 0 are in fact highly asymmetric tops {Ib and 
Ic differ by a factor of 3). In the absence of any 
better procedure, we assume that H 2 0 and D 2 0 can 
be treated as slightly asymmetric tops and find that 
they must then be classified as near-spherical tops. 
For such molecules the resonant correction is 

5=—/ )( 3 gjv* p (wn ~ 7/2 

16 \16 a ok t) Ci n t I a (IbIc) *1.4 4 I A J ' 

(7) 

[Note that Eq. (7) reduces to Eq. (3) if the moments 
of inertia are equal.] From Eq. (7) and Table IV (and 


Table VI. Thermal conductivities for analysis of the data on 
H 2 0-D 2 0 mixtures (/real cm -1 sec -1 °K“ 1 )* 



Temperature, °K 

Gas 

478.0 

525.6 

H s O, (exptJ.) 

79.5 

89.4 

H 2 0, nonresonant (calc.) 

98.3 

111.1 

D 2 0, (exptl.) 

80.0 

90.2 

D 2 0, nonresonant (calc.) 

88.8 

100.3 

Xj 2 , nonresonant [calc., Eq. (11)] 

93.7 

105.8 

Equimolar H 2 0-D 2 0 [calc., Eq. (10)] 

86.2 

97.1 

Equimolar H 2 0-D 2 0 (exptl.) 

79.8 

89.9 


assuming the diffusion collision cross section for D 2 0 
is also 3% larger than for H 2 0) we find 

5d 2 o/5h 2 0 = 0.3613Cmt H2 o/ C intx) 2 0* (8) 

Equations (1) and (4) suggest a simple average 
moment of inertia may be used to calculate a mean 
ZrotT 1 . Hence we have assumed 

Zd 2 o/Zh 2 o = 0.5933. (9) 

Finally, the dimensionless quantity aD n /v involv- 
ing the self-diffusion coefficient was computed as 
(6/5) ( A *). 22 The quantity (^4*) is a ratio of collision 
integrals (essentially the ratio of the viscosity and 
diffusion cross sections) and is quite insensitive to 
temperature and the details of the intermolecular force 
law, at least for spherically symmetric potentials. In 
the absence of experimental information on self-diffu- 
sion in water vapor, we have taken (A * ) from calcu- 
lations for a modified (angle-independent) Stockmayer 
potential. 23 This is the only place in the analysis where 
we have had to make any assumption about the inter- 
molecular force law. 

Collision numbers for rotational relaxation and reso- 
nant corrections deduced from the data analysis are 
presented in Table V. The collision numbers for H 2 0 
(one to three collisions) seem low for a molecule with 
such small moments of inertia. (Methane, for example, 
requires perhaps 10-15 collisions. 8 - 9 ) However, these 
low values can perhaps be rationalized by noting that 
the large dipole moment causes a large deviation from 
spherical symmetry in the intermolecular force field 

22 Reference 10, Eq. (8.2-48), p. 540. 

23 L. Monchick and E. A. Mason, J. Chem. Phys. 35, 1676 
(1961). 
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and consequently promotes energy transfer between 
rotation and translation. 

The resonant corrections are seen to be somewhat 
smaller than values computed from Eq. (7), especially 
at the lower temperatures. The derivation of Eq. (7) 
involves among other approximations the replacement 
of summations by integrations. In the case of ammo- 
nia, Monchick and Mason 1 carried out the direct sum- 
mation and found that the analytic expression for the 
resonant correction gave too high a value at low tem- 
perature. A similar deviation for water might have 
been anticipated. All in all, the data on pure H 2 0 and 
D 2 0 tend to confirm Mason and Monchick’s postulate 
regarding the importance of resonant transfer of rota- 
tional quanta in reducing the heat conductivity of 
highly polar gases. 

Unfortunately, the data on the H 2 0-D 2 0 mixtures 
seem to negate this conclusion. It is a simple matter 
to compute the thermal conductivities of H 2 0 and D 2 0 
using the collision numbers of Table V but omitting 
the resonant correction (values so computed are shown 
in Table VI). The effect of the resonant correction is 
substantial since it reduces the calculated conductivity 
of D 2 0 by 10% and H 2 0 by almost 20%. 

In H 2 0-D 2 0 mixtures the collisions between H 2 0 
and D 2 0 molecules should not involve resonant ex- 
change and should consequently be characterized by 
these nonresonant conductivities. Although we have 
no rigorous mixture theory for polyatomic gases, we 
may make a reasonable prediction based on monatomic 
gas theory. The rigorous expression 24 for a binary mix- 
ture of isotopic monatomic gases of nearly equal mass 
is, to a very good approximation, 


Xi Xi 2 X 2 

Here xi and x 2 are the mole fractions of the component 
gases, while Xi and X 2 are their thermal conductivities. 
The quantity X i2 is defined by Eq. (8.2-34) of Ref. 10. 
For isotopic species the equation reduces to 

Xi2 = Ci (Xi 2 +X 2 2 )J*. (11) 

Clearly, the first and third terms of Eq. (10) corre- 
spond to self-collisions, whereas the second term arises 
from collisions between the unlike species. Thus it 
would seem a reasonable approach to compute X i2 from 
Eq. (11) using the nonresonant conductivities of Table 
VI. Thermal conductivities of the equimolar mixtures 
can then be obtained from Eq. ( 10) using the measured 
values for H 2 0 and D 2 0. Results in Table VI show the 
computed mixture conductivities exceed the experi- 
mental values by about 8%. This is far beyond the 
range of any conceivable experimental uncertainty. 

24 Reference 10, Eq. (8.2-36), p. 535. 


As a matter of fact, the experimental mixture con- 
ductivities fall almost exactly midway between the 
conductivities of the pure components. 

In our opinion, this suggests rather strongly that a 
resonant exchange of rotational quanta does not play 
an important role in determining the thermal conduc- 
tivity of water vapor. Admittedly Eq. (10) is not 
rigorous for mixtures of polyatomic gases (a correct 
formulation would be considerably more complicated 
and involve, among other things, relaxation times 
among the like and unlike molecules explicitly). None- 
theless, we feel that Eq. (10) must be a good ap- 
proximation for mixtures of isotopically substituted 
polyatomic gases because of its simple collisional in- 
terpretation. In computing \ 12 from fictitous nonreso- 
nant conductivities of the pure gases, we are tacitly 
assuming that the collision number for D 2 0 relaxing 
H 2 0 is the same as for the self-relaxation of H 2 0 (with 
a similar assumption for H 2 0 relaxing D 2 0). This is 
probably not strictly correct, but if Eq. (4) is any 
guide, the assumption is in error by no more than 
10% (due to the mass difference — more likely 5%, 
since an average mass is probably appropriate). It 
would seem to us to be a most remarkable coincidence 
indeed if the collision numbers for the unlike interac- 
tions were smaller than those for self-collisions so as 
to exactly compensate for absence of resonant exchange 
of rotational quanta in H 2 0-D 2 0 collisions. • 

CONCLUDING REMARKS 

Although the experimental conductivities of pure 
H 2 0 and D 2 0 are consistent with the postulate 1 that 
a resonant exchange of rotational quanta is important 
in determining the thermal conductivity of polar gases, 
the data on the equimolar H 2 0-D 2 0 mixture seem to 
refute this idea. Nonetheless, the fact that the con- 
ductivity of water vapor appears anomalously low in 
relation to its viscosity needs an explanation. Tenta- 
tively, we suggest that exchange of rotational energy 
on grazing collisions is indeed important, but it is 
essentially classical rather than quantum mechanical. 

Since a large dipole moment causes a long-range asym- 
metry in the intermolecular force field, it should pro- 
mote such exchange. However, if the phenomenon is 
classical, exact matching of rotational energy levels is 
not a requirement, and exchange between dissimilar 
polar molecules should occur. 

These conclusions are not definitive because of two 
theoretical shortcomings: we do not as yet have ade- 
quate theories for asymmetric top molecules or polar 
gas mixtures. The first shortcoming is overcome most 
easily by studying gas pairs for which the theory does 
exist: linear dipoles, such as HC1 and DC1, and symmet- 
ric tops, such as NH 3 and ND 3 . We hope to initiate 
such experiments shortly. In the meantime, perhaps a 
theory for gas mixtures will be forthcoming. 



